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Abstract 

Many models possess unwanted relics, which should be diluted by 
entropy production just before the big-bang nucleosynthesis. A field 
responsible for the entropy production may produce stable weakly in- 
teracting massive particles, if kinematically accessible. We compute 
their relic abundances by integrating out coupled equations numer- 
ically. Applying our results to supersymmetric standard models, we 
argue that the neutralino lightest superparticles will overdose the uni- 
verse in most of the parameter space if the reheat temperature is of 
the order of 10 MeV. 



Many models beyond the standard model predict unwanted relics, which 
should be diluted by entropy production such as inflation and subsequent 
reheating. Among other things, weakly interacting massive scalar fields will 
be the most problematic, as they can store their energy density in the form 
of coherent oscillation, which may not be diluted by inflation because the 
coherent mode during the inflationary epoch is generally displaced from its 
value in the true vacuum. An example is a problem associated with moduli 
fields of string theory, named the cosmological moduli problem In string 
theory coupling constants are determined dynamically by vacuum expecta- 
tion values of the moduli (and dilaton) fields. These fields generically acquire 
masses comparable to the gravitino mass, which should be in the TeV range 
in gravity mediated supersymmetry breaking. Damped coherent oscillation 
and subsequent decay of the fields would spoil the standard big-bang nu- 
cleosynthesis (BBN). Late-time entropy production before the BBN should 
operate to dilute the energy density of the problematic scalar fields. Indeed 
thermal inflation was proposed 0, [J, for which reheat temperature should 
be much lower than the electroweak scale for a successful dilution. 

Another example of the problematic scalar fields is an invisible axion with 
decay constant larger than 10 13 GeV, whose coherent oscillation would exceed 
the critical density of the universe ||. Entropy production may be able to 
dilute the energy density of the axion field ||. Since the oscillation of the 
axion field starts around the QCD phase transition, the reheat temperature 
at the entropy production should be lower than 100 MeV. Interestingly it 
has been argued that strongly coupled Eg x Eg heterotic string theory (M- 
theory) possesses such an axion candidate with the decay constant ~ 10 16 
GeV, provided that non-perturbative effects to generate the superpotential 
for the axion field are suppressed [|, |7|. To dilute the energy density of the 
M-theory axion to a harmless level, the reheat temperature of the entropy 
production should be at its lowest value allowed by the BBN, 1-10 MeV. 

The requisite entropy production is provided by decays of massive unsta- 
ble particles. The decays, however, may produce too many stable particles 
as well if the production is kinematically allowed. In a supersymmetric the- 
ory, the lightest superparticle (LSP), which is likely to be a neutralino, is 
stable under the assumption of i?-parity conservation. Indeed it was argued 
that the neutralinos produced by the decays of the massive particles tend to 
overdose the universe, if the reheat temperature is much lower than 1 GeV 

Bl- 

Motivated by the aforementioned demands for the late-time entropy pro- 
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duction mechanism, we shall reexamine the problem of the relic abundances 
of the neutralinos which are produced by the decays of massive unstable par- 
ticles. Given an annihilation cross section of the neutralinos, we will solve 
coupled equations numerically to evaluate the relic abundances, taking into 
account effects of exponential decays of the unstable particles, whereas an 
approximation of simultaneous decay was made in the previous works || |9| . 
It turns out that results obtained in our numerical calculation agree with 
the previous estimates up to a numerical factor of order unity and hence we 
confirm the validity of arguments in Refs. [§, |J . Note that our computation 
should be applicable to any stable weakly-interacting-massive-particles. We 
will then argue that the stable LSP suffers from the overclosure problem in 
most of the parameter space if the reheat temperature is of the order of 10 
MeV or less. We will also briefly mention possible implications of our anal- 
ysis when there remains the M-theory axion in the massless spectrum of the 
heterotic M-theory. 

We consider the situation that massive unstable particles (f) of the decay 
width T come to dominate the energy density of the universe, then decay 
to neutralinos \ as well as to radiation with branching ratios B and 1 — B, 
respectively. By assuming two-body decay, the energy density of x j us t 
after the production is They subsequently lose their energy, by the 

scattering with thermal bath, down to ~ m x T). Assuming this energy 
transportation occur rapidly, we define the effective branching ratio B e ^ = 

£> m X 

Time evolution of cosmic scale factor R and energy densities p^, pn, p x of 
0, radiation, and Xi respectively, are described by the following Friedmann 
equation and (integrated) Boltzmann equations: 

+ PR + Px) > (!) 

(2) 

dt — 3 d t ^-^)r(^) + ^(^) 2 ^, 0) 

^p x ) = B cS T(R3p,) - ^Vpx) 2 ^ , (4) 

where M is the reduced Planck mass: it is related to the Newton's constant 
as G = 1/wip! = l/(87rM 2 ), and (av) is an ensemble-averaged annihilation 
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cross section of neutralinos times their relative velocity. We treat (erf) as a 
constant, which is indeed the case when the S"-wave is dominant. When a 
higher partial wave dominates, on the other hand, (av) should be interpreted 
as a sort of weighted average in the time interval around when the annihi- 
lation process is important. In the above equations, we did not include pair 
production of the neutralinos, which is not effective at low temperatures. 
Note that the right-hand side of the Boltzmann equation for p R is consistent 
with the energy conservation within comoving volume for the case R = 0. 
Following [|Kj we introduce the dimensionless variables: 

x = Tt , z = R/Ri , 

U = P<t>/ P<t>i J fR = PR/ P<t>i > fx = Pxl P<t>i J 



F^z'U, F R = z 4 f R , F x = z 3 f : 



x H = THr 1 = T 



m x T 
(av) 



3M 2 
fxo = Pxo/ P<t> 



l + fm + f- 



X ' 

1-1/2 



X'/ 



(5) 
(6) 
(7) 

(8) 
(9) 



Here subscripts i mean their initial values. By using these variables, the 
coupled differential equations above can be written as follows: 



z 
z 



x H (1 + fm + f xi 

F'j, = —F^ , 
1 



-1/2, 



z~ 3 F + z~ 4 F R + z~ 3 F x 



,1/2 



F' R 
F 'x 



B eS )zF^ + f^z~ 2 Fl 



BcsF^ 



f-l -3rf 



(10) 

(11) 

(12) 
(13) 



'x0~ -x 

where prime denotes d/dx. We impose the following initial conditions at 

X Xi , 



1 



Pm 



F Rl = — < 1 



= 1 , 

F xi = (« 1) 



(14) 
(15) 



There is no difficulty in solving the Boltzmann equation for 

F^ = e-C*-**) . 
By assuming p oc R~ 3 for < x < 1, 

and TH- 1 



y 3/2 _ ^ X 



2x H 



3 



(16) 



(17) 
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for Xj <C x < 1. Our numerical calculation shows that z\ x= i is about 10% 
larger than the estimate (t^ff 1 ) 2 ■ 

With an approximation of the simultaneous <fi decay at x — 1, one obtains 

2/3 



2 *H 



1.20 x (1.09 x^ 2/3 ) . (18) 



Numerical calculation |10| for the coupled equations shows, on the other 
hand, 

roo 

Fr\x>i = / ze' x dx 
Jo 

= l.09x H 2/3 , (19) 

for Fjh < 10 -3 and B c r < 10 -3 . We confirm that the formula above is correct 
also for larger Fm and/or B e R with at most about 10% error. 

Our main concern in this paper is F x for x ^> 1, to which p x is propor- 
tional. A crude behavior of F x can be seen as follows. If f~ l 
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negligibly small, or quantitatively if B eS x 2 H f x ^ < 1 as we shall see shortly, 

_ J B eS (x - Xi) + F xi , for Xi < x < 1 ( . 

x ~\B cS , for x > 1 . 

That is, F x increases monotonically with x. On the other hand, if is 
not so small the second term of the Boltzmann equation for F x cannot be 
neglected when F x becomes large: 



F' 

x 



B^e^l - 1.2B eS x 2 H fJ] . (21) 



Thus when B eS x 2 H f x ^ > 1, F x increases until x ~ 1, but decreases thereafter 
to reach F x \ x:> i < B eS . Note that 

Besxy-J = 6M 2 BT^- 

= 117 (lO-iOGeV- 2 ) (lb^Gev) (as) (te) ' (22) 
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Fig. 1 shows the evolution of F x , obtained by integrating out the coupled 
differential equations numerically. It shows it increases till x = 2 or so, then 
decreases to reach a constant value. In the same figure, we also show the 
evolutions of the other quantities, z, F$ and Fr. 

Let us now compare our numerical results with a naive estimate given 
in |J. This naive estimation of i ? x | a; >i for the case where the annihilation 
process becomes effective was obtained by equating the annihilation rate 
(av) p x / m x t° the expansion rate H at the decay epoch (x = 1) of 0, with 
the approximation of the simultaneous decay. One finds 



m x H 

Px\x=l 



(av) 



x=l 



(23) 

3 {av) ' { } 



or, 

F x \ x ^>i = F x \ x= i = -x H f x0 . (24) 

On the other hand, we find the following fit for our numerical results (see 
Fig. 2): 

F x \ x>>1 = 1.52(B cS x 2 H f^y°- 9W B cS (25) 

for 10 < B cS x 2 H f x J < 10 8 within 15% error. Comparison of the two results 
shows that the naive estimate is correct up to a factor of order unity in a 
wide range of the parameters. As B eS x 2 H f~Q oc BT^j- becomes larger, F x \ x:> i 
is more underestimated if one uses the naive estimation. For example, for 
B cS = 1CT 1 and x 2 H f x ^ = 10 8 , 

FX^ = 1.5 x (10 8 )- 1 = 1.5 x 1(T 8 , (26) 



while, 



^numerical = ^ x ^ Q -! x ^8^0.919 x 1Q -1 



= 5.61 x 1(T 8 

= 3.7xF x C™. (27) 
Now we turn to p x /s\ x - > i where s stands for entropy density. Using 

p R = aT 4 , 

8 = % = -a^pT , (28) 
3 T 3 HR ' K ' 

2 

a= 30 9 * 
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with g* being the effective relativistic degrees of freedom, we find 



Px _ 3 P<l>i F X 
s 4 a V*F% 4 



3 5/4 



4 01/4(^3^)3/4 



(A4T) 



1/2 



(29) 



Substituting our numerical results Eqs. (pT9|) and (|25|) , we obtain the following 
fit 



Px 
s 



35/4 



x>l 



1.09) 



-3/4 



7T 



-1/4 



30* 



xl.52 6M 2 fir 



3.08/T 2 x 10 4 P — 



(avY 



(Mr) 

-0.919 



1/2 



(av) 



-0.919 



X 



Vio- lo Gev- 2 ; 

/ m \ -°- 081 / B 



V10 3 GeV 



V10 2 GeV 

0.081 



_l \ 0.419 



0.5 



9* 



lsec J 

0.46 



10.75 



(30) 



where p C r/s|o = 3.64/i 2 x 10~ 9 GeV is the ratio of the critical density p cr to 
the entropy density s at present. This formula reproduces values of p x /s| x> i 
calculated numerically within 15% error for 10 < B ca x 2 H f-J < 10 8 . Alterna- 
tively, the annihilation cross section is related to the relic abundance as 



(av) = 0.767/r 2 ' 18 x lO^GeV -2 

-1.09 



X 



Pxl s \x 



>1 



Pcr/<s|o 



10 2 GeV 



1.09 



-i-l 



lsec 



0.456 



x 



/ ///,.. \ °- 088] / n • 



V10 3 GeV/ 



0.5 



9* 



10.75 



-0.501 



One finds that Eq. (31) agrees with the previous estimate in Ref. 
a factor of order unity, if one recalls the temperature when <p decays is 

1 \ -1/2 



-A/2 

T ~ g- 1/A ^MT ~ MeV 

\ lsec / 



(31) 
up to 

(32) 



which may be referred to as the reheat temperature. 
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The ensemble- averaged annihilation cross section (times relative velocity) 
(av) can be, in general, expanded in the averaged velocity squared (v 2 ) of 
the LSP: 

(av) = a + b(v 2 ) + ■■■ . (33) 
The coefficients a, b, etc. strongly depend on the composition of the LSP as 



well as superparticle mass spectra [11 . Sometimes the coefficient a of the 



S'-wave becomes very small due to the Majorana nature of the neutralino. 
This is indeed the case when the LSP is the bino-dominant neutralino, in 
which a is proportional to the mass squared of the fermion in the final state. 
Then, in order to evaluate (av) and subsequently Q x h 2 , we have to know 
(v 2 ) precisely. The authors of Ref. || found that the bino-dominant LSP is 
likely to reach kinetic equilibrium and thus (v 2 ) ~ 3T/m x when the reheat 
temperature T is somewhat high, T > a few x 10 MeV. Using this value of 
(v 2 ), they concluded that the LSP would overdose the universe as far as the 
reheat temperature is lower than a few hundred MeV. It is not clear, however, 
whether the LSP is in the equilibrium or not for other sets of parameters. 
Without detailed knowledge of the velocity distribution, there are still some 
conclusions we can draw in general. For instance, in view of the fact that 
the annihilation cross section is generally bounded from above as 

1CT 2 . . 

M £ — ^ T 34 

X A 

we may conclude that the LSP with the mass at a few tens GeV or more 
does not annihilate sufficiently enough to survive the overclosure limit, if the 
reheat temperature is relatively low, say T < O(10) MeV, or T" 1 > (9(10~ 4 ) 
sec, which seems to be required to make the M-theory axion viable. Note 
that the bound does not apply when s-channel resonances (by Z boson 
or Higgs bosons) occur. To discuss the relic abundance in this case, it seems 
that we need to evaluate the velocity distribution of the LSP, which is beyond 
the scope of the present paper. 

Finally we would like to discuss implications of our analysis to the M- 
theory if it turns out that the axion with large decay constant indeed arises. 
One needs to invoke an entropy production mechanism to dilute the energy 
density of the axion's coherent oscillation, unless its initial amplitude can be 
tuned to be very small. As we argued, the entropy production would produce 
too many neutralinos, if kinematically allowed, which would not annihilate 
efficiently. As a result they would overdose the universe, if they remain today. 



7 



This suggests that R-parity should be violated, making the neutralino-LSP 
unstable. Of course there are many loopholes to this argument, and we 
would not insist that the R-parity must be broken in the heterotic M-theory. 
Rather we regard the R-parity violation as in interesting option to which 
more attention should be paid. 
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Figure captions 

Figure 1 Time evolution (x = Ft) of z = R/Ri, F^ = z 3 f^, F R = z A f R , 
F x = -2 3 / x , and p x / s, where /* = p*/ p^j (* = 0, it!, x) and T is the decay 
width of 0, for the case when B c r = 2m x /m ( p = 10 _1 , xh = FHf 1 = 
10" 2 , / x0 = m x T/((av)p < j )0 ) = 10~ 9 , T" 1 = lsec, and g* = 10.75. We 
can see that entropy S oc F^ 4 within comoving volume increases during 
0.01 < x < 2. We can also see F x increases until x ~ 2 but decreases 
thereafter to reach certain constant value. Likewise, p x / s decreases for 
x > 2 to reach certain asymptotic value. 

Figure 2 Asymptotic values of F x and p x /s as a function of f~Q, for the 
case _B c fj = 10 _1 , xh = 10~ 2 , = lsec, and g* = 10.75. Qualitative 
behavior of F x \ x ^i and p x /s|x>i changes depending on if B e QX 2 H f~Q is 
greater than 1. We also show naive estimations, -F x |°>\ e = (3/2)x^ 2 / x o 
and p x /s|£»T = 2~ 3 / 2 3 3 / 4 a- 1 / 4 v / MT ^ X |£>>T, and linear fits for 10 < 
B c ftX 2 H f~Q < 10 8 , for comparison. Linear fits are accurate within 15% 
error. 
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